Structure and stability of a high-coverage (1x1) oxygen phase on Ru (0001) 
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The formation of chemisorbed 0-phases on Ru (0001) by exposure to O2 at low pressures is apparently limited to 
coverages O < 0.5. Using low-energy electron diffraction and density functional theory we show that this restriction 
is caused by kinetic hindering and that a dense O overlayer (0 = 1) can be formed with a (1x1) periodicity. The 
structural and energetic properties of this new adsorbate phase are analyzed and discussed in view of attempts 
to bridge the so-called "pressure gap" in heterogeneous catalysis. It is argued that the identified system actuates 
the unusually high rate of oxidizing reactions at Ru surfaces under high oxygen pressure conditions. 

PACS numbers: 68.35.Bs, 61.14.Hg 



Chemisorption studies with well-defined single crystal surfaces are propelled by the prospect to gain deeper insight 
into the elementary steps governing heterogeneously catalyzed reactions. Experiments of this kind are usually con- 
ducted under ultrahigh vacuum (UHV) conditions with partial pressures typically < 10~^ mbar in order to apply 
the standard tools of surface physics and to control the state of the surface. In contrast, "real" catalytic reactions 

■ take place at much higher pressures (frequently even exceeding 1 bar); this difference is sometimes referred to as the 
"pressure gap" . 

Although there exists a number of examples for which extrapolation of data over a wide pressure range has been 
^ I demonstrated to be safely justified Q, such a conclusion may by no means be generalized. The concentrations 

■ of certain surface species relevant for catalysis at high pressures may be so small under low pressure conditions 
that they escape detection, either for thermodynamic (low stability and hence short surface residence time) or for 
kinetic reasons (low sticking coefficient). Such an example is presumably the Ru(OOOl) surface interacting with O2: 

' Dissociative chemisorption under low-pressure conditions leads to apparent saturation at a coverage Q = 0.5 (the 
\^ ■ coverage Q being defined as the ratio of the concentration of adparticles to that of substrate atoms in the topmost 
0^ ' layer) associated with the formation of an ordered (2x l)-superstructure Since the smallest separation of O 

atoms in the (2xl)-0 overlayer is only one lattice constant, also a (Ixl)-O phase with 8 = 1 in which the O atoms 
form the same lattice as the Ru(OOOl) seems to be feasible (at least from purely geometrical arguments). Indeed, 
from high pressure experiments of the catalytic oxidation of carbon monoxide on Ru (0001) the formation of a high- 
coverage 0-phase was speculated Q , and recent density-functional theory (DFT) calculations suggested the existence 
of an ordered O adlayer with 8 = 1, for which the oxygen is still well chemisorbed, but with a much weaker bond 
than for the other, more open phases, making the 8 = 1 adlayer particularly effective for oxidizing reactions An 
experimental proof of this theoretical result, and in particular, an analysis of whether or not oxygen might be present 
below the surface Ru layer is still missing. If the formation of the high coverage phase were indeed only kinetically 
^ . obstructed, then it may be realized by high oxygen pressure or by exposure to a more strongly oxidizing molecule such 
■ ■ ' as NO2. The latter approach was first employed to reach high O coverages on Pt (111) ^,0, and more recently it was 
demonstrated that also for Ru(OOOl) a concentration of surface oxygen can be obtained reaching, or even exceeding, 
8 = 1 |,|. 

In the present letter we report the confirmation of the formation of a high-coverage 8 = 1 oxygen phase on 
Ru (0001) where, through the methods of low-energy electron diffraction (LEED) and density-functional theory (DFT) 
calculations, we determined the detailed atomic geometry and the associated energetics which give insight into the 
behavior of O on Ru (0001) in general. Using the technique of dissociative NO2 chemisorption, we show that this 
phase is indeed an on-surface (1x1) adlayer with a negligible amount of subsurface oxygen. 

The experiments were conducted in a UHV chamber (base pressure = 1 x 10^^*^ mbar) equipped with a display 
type four grid LEED optics and with standard facilities for surface cleaning and characterization. Details about the 
experimental set-up and sample preparation can be found elsewhere pO| ]. The preparation of the oxygen overlayer 
with 8 = 1 was accomplished by first exposing the Ru(OOOl) surface to O2 at room temperature until saturation 
was reached, i.e. the (2xl)-0 phase was formed. Subsequently, the Ru (0001)-(2x l)-0 surface was exposed for 15 
minutes to NO2 at a pressure of 5 x 10^^ mbar. The sample temperature was chosen now to be 600 K so that during 
decomposition of NO2 into adsorbed O and NO (the dissociative sticking coefhcient of NO2 is about one), only O 
remained on the surface while NO was instantly released into the gas phase, thus enriching the oxygen content at the 
surface; note that NO desorption from Ru (0001) is completed at about 500 K ||ll|. Concomitant LEED measurements 
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indicated at first a sharpening of the (2x2) pattern followed by a gradual transformation into a bright (Ixl)-pattern 
with low back-ground intensity, consistent with the formation of a (Ixl)-O phase. This (Ixl)-O structure remained 
stable even upon further NO2 exposure so that no additional oxygen accumulated at the surface, either because the 
sticking coefScient dropped or the excess oxygen atoms disappear into the subsurface region. Subsequently recorded 
thermal desorption spectra of O, in combination with the structural analysis (described below), give compelling 
evidence that additional O atoms are indeed formed which dissolve into the bulk of the sample, leaving the oxygen 
concentration in the (Ixl)-O overlayer constant. A similar behavior was also observed by Weinberg LEED 
intensity measurements were performed at normal incidence of the primary beam at a sample temperature of 100 K. 
A computer-controlled video camera was used to record integrated spot intensities of five integral-order beams (energy 
range 50 to 620 eV) from the fluorescence screen. LEED curves were calculated using the program code of Moritz JlJ] 
and compared with the experimental LEED curves by applying an automated least-squares optimization scheme [nSf, 
based either on the reliability (i?) factors i?DE [0 or Rp p^ ]. 

For the DFT calculations we employed the generalized gradient approximation (GGA) of Perdew et al. for 
the exchange-correlation functional and used ab initio, fully separable pseudopotentials. The surface was modeled by 
a supercell consisting of four layers of Ru(OOOl) where the O atoms were adsorbed on one side of this slab |]l7| , |l8| . 
We relaxed the positions of the O atoms and of the Ru atoms in the top two layers, keeping the lower two Ru layer 
spacings fixed at the bulk values. Details of the calculations, as well as a discussion of the clean Ru surface and more 
open O adlayers, can be found in Ref. js). 

In the LEED analysis, the 0-Ru layer spacing and the first two Ru-Ru layer spacings of the Ru (OOOl)-(l x l)-0 
surface were simultaneously refined, starting the automated search from different adsorption sites of oxygen, namely 
the fee-, hep-, on-top- and bridge-sites. The resulting optimum i?- factors for the quoted adsorption sites are compiled 
in Tab. I from which it becomes evident that the hep-position is clearly favored. In addition to these standard con- 
figurations, we tested a (1x1) structure placing O in both the fee and hep position with variable concentrations. The 
mixing of adsorption sites with variable concentrations was simulated within the framework of the averaged-t-matrix 
approximation | |T9[ |. The multiple scattering between neighboring 0-atoms in fee and hep-sites could alternatively be 
switched off. The best agreement with the experiment was in all cases achieved by the adsorption of O atoms exclu- 
sively in hep-sites. Even a small concentration (0 — 0.1) of 0-atoms in fee sites deteriorated the theory-experiment 
agreement appreciably, as indicated by an increase of the i?- factors by 0.05. The next issue we addressed carefully 
was as to whether the results are sensitive to the presence of sub-surface oxygen. For this purpose we put, in ad- 
dition to the O-(lxl) overlayer, O atoms, with variable concentrations, into the octahedral sites between the first 
two Ru-layers and optimized their concentration and the structural parameters. The optimization scheme always 
ended with a structure in which O adsorbs again solely in hep-sites on the surface. On the basis of these simulations 
we can therefore safely rule out any O concentration exceeding 0.1 monolayer (ML) between the top two Ru layers. 
Thermal desorption spectra (TDS) of oxygen, on the other hand, indicate that there is much more oxygen than = 1 
present in this system and that appreciable amounts of oxygen are still present on the surface even after many thermal 
desorption cycles. Obviously, these O atoms prefer to dissolve into the bulk rather than to form a subsurface phase 
between the first two Ru-layers. 

The optimum atomic geometry of the Ru (0001)-(lxl)-O phase, as provided by the LEED analysis, is displayed in 
Fig. 0. The intensity spectra for this best-fit geometry are shown in Fig. ||; the corresponding i?- factors are i?DE = 0.21 
and Rp — 0.23. The Ru-0 bond length of 2.00 A ±0.03 A is similar to the data reported previously for the (2x2)-0 
(0 = 0.25) and the (2xl)-0 (0 = 0.5) structures which are 2.03 A and 2.02 A respectively [|o). A remarkable 
feature of the (Ixl)-O phase is the substantial expansion of the topmost Ru layer spacing (2.21 A) by 3.7 %. This 
expansion is comparable to that derived for the O/Zr(0001) system (3.1 %) where IML of oxygen is reported to be 
evenly distributed in octahedral sites between the first and second, and second and third Zr layers pl[| . 

The experimentally derived structural data agree nicely with the results of total energy calculations performed for 
varying 0-coverages, namely the (2x2)-0, (2xl)-0, and (Ixl)-O structures. For all these phases, the energetically 
most favorable adsorption site was found to be the hep-site. In Tab. II, the structural parameters for the (Ixl)-O 
system provided by DFT-GGA and LEED are compared; the agreement convincingly demonstrates the power and 
reliability of the DFT-GGA calculations. It is important to notice that the large expansion of the top Ru interlayer 
spacing is not due to subsurface oxygen. It results from the binding of the surface Ru atoms with on-surface oxygen 
which empties bonding Ru-Ru d states, weakening the attraction between the top and second substrate layers. In 
Fig. ^, the calculated binding energies for an oxygen atom in these three oxygen phases are compiled for the cases 
where the O atoms reside in hep and in fee positions. It can be seen that with increasing 0-coverage the binding energy 
becomes smaller which refiects a repulsive interaction between the adsorbates and is stronger for oxygen in the hcp- 
than in the fee-site. Using these results for a simulation of thermal desorption spectra we obtain good agreement 
with experimental data corresponding to the coverage range < 0.5. For oxygen in the hep-site at = 1, the 



2 



energy required to remove one oxygen atom, i.e. the energy to create an oxygen vacancy in the adlayer is quite small, 
namely, 1.2 eV. This result was estimated from calculations with a (2 x 2) surface cell (the artificial vacancy-vacancy 
interaction is negligible for our concerns, i.e., less than ±0.15 eV, as tested by larger surface cells). Comparing the 
binding energies per atom of the two = 1 layers, we find that the energy difference is only 0.06 eV. Because the 
= 1 layer forms subsequently after the completion of the lower coverage phases for which the hep site is by several 
tenths of an eV more favorable, it has been argued ||^ that, despite the small energy difference, the = 1 phase might 
be a rather perfect hep-site adlayer. Nevertheless, it is obvious that an independent, experimental structure analysis 
of this system was mandatory. The very good agreement between the structural parameters obtained by DFT-GGA 
and LEED, discussed above, provides the required confirmation. 

The above results imply that the saturation 0-coverage Q ~ 0.5, which is reached after exposing a Ru (0001) surface 
to O2 at low pressure, is only apparent and in fact limited by the kinetics of dissociative adsorption. Since the rate of 
impingement increases in proportion to the partial pressure, catalysis under high pressure conditions with an excess 
of O2 in the feed becomes rather likely to involve surface phases with 0-coverages approaching 8 = 1. Given the fact 
that the metal-oxygen bond strength is rather weak at the full coverage (see Fig. ^), it is expected to improve the 
reactivity for reactions of the type Oad + X — > OX. The present system hence represents, in our opinion, an example 
for which extrapolation of the kinetics of a catalytic reaction across the "pressure gap" may be questionable, as also 
suggested by the quoted work on CO oxidation Q . 

Extrapolation of the data of Fig. || to even higher coverages suggests that adsorption of oxygen even beyond 8 = 1 
might still be exothermic. DFT-GGA calculations at coverages 8 = 1.25 indicate, however, that further incorporation 
of O chemisorbed on the surface into the (Ixl)-O structure is unstable. On the other hand, the occupation of the 
subsurface octahedral adsorption site, just below the first Ru layer, is still exothermic; it is however, appreciably less 
favorable than on-surface sites for coverages 8 < 1. Interestingly, for even deeper positions, the binding energy of 
these additional O atoms is similar. Furthermore, the diffusion energy barriers for oxygen penetrating through the 
first and second Ru layers are found to be such that the barrier for passing through the second layer is significantly 
smaller than that of the first, O-covered Ru layer. As a consequence, the theoretical results suggest that subsurface 
oxygen will dissolve into the bulk on entropy grounds. ^From these calculations (and also from our TDS experiments) 
the existence of a close-packed O overlayer with 8 = 3 on Ru (0001) as proposed by Malik and Hrbek ^ can clearly 
be ruled out. Hence, an impinging NO2 molecule gets cither reflected by the (Ixl)-O surface or dissociation takes 
place and the oxygen atom is incorporated into the subsurface region from where it proceeds into the bulk. 

In summary, from a combined investigation using DFT-GGA and LEED we confirmed the existence of a high- 
coverage (8 = 1) phase of oxygen on Ru(OOOl) and provide a detailed identification of the geometry. Furthermore 
we provided evidence that atomic oxygen may enter Ru-bulk once all the on-surface hep sites are occupied, but not 
staying between the first and second substrate layer. This unusually high coverage oxygen structure on Ru (0001) 
is likely to enhance the rate of oxidation reactions (e.g. of carbon monoxide) and other surface reactions involving 
adsorbed oxygen. It also immediately raises the questions of the possible existence of other high coverage surface 
structures likewise achievable by bypassing the "pressure gap" . 
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O-adsorption site 


Rp 


Rue 


on top 


0.68 


0.47 


bridge 


0.60 


0.39 


fee 


0.59 


0.49 


hep 


0.23 


0.21 


no oxygen 


0.55 


0.37 



TABLE I. Optimum i?-factors for different structural models of Ru(0001)-(lxl)-O. 



parameter 


DFT-GGA 


LEED 


Ru-O layer spacing (A) 


1.26 


1.25 ± 0.02 


Ru-O bond length (A) 


2.03 


2.00 ± 0.03 


first Ru interlayer spacing (%) 


+ 2.7 


+ 3.7 ± 1.4 


second Ru interlayer spacing (%) 


-0.9 


-0.5 ± 1.8 



TABLE II. Structural parameters of the (Ixl)-O phase on Ru(OOOl) as obtained by DFT-GGA calculations (see text and 
ref. 4) and the LEED analysis. 



4 



FIG. 1. Top view (a) and side view (b) of the atomic geometry of (lxl)-0/Ru (0001) with O sitting in the hcp-hollow site. 
The arrows indicate the direction of the displacements of the substrate atoms with respect to the bulk terminated positions. 
The distances are those determined by LEED and are given in A. 
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FIG. 2. Comparison between experimental and calculated LEED curves for the best-fit geometry of the (Ixl)-O phase on 
Ru(OOOl). The overall i?-factors are Rue = 0.21 and Rp = 0.23. 

FIG. 2 STAMPFL ET AL. 
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FIG. 3. Binding energy of O on Ru(OOOl) at various coverages, with respect to 1/2 O2 molecule. Adlayers are assumed with 
fcc-hoUow adsorption sites (dashed line) and with hep- hollow sites (continuous line). 
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